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Abstract 

We study a bottom-up holographic model of large- N c Yang-Mills theory, in which conformal 
invariance is broken through the introduction of a dilaton potential on the gravity side. We use 
the model to calculate the spectral densities of the shear and bulk channels at finite temperature. 
In the shear channel, we compare our results to those derived in strongly coupled Af = 4 Super 
Yang-Mills theory as well as in weakly coupled ordinary Yang-Mills. In the bulk channel, where 
the conformal result is trivial, we make comparisons with both perturbative and lattice QCD. 
In both channels, we pay particular attention into the effects of conformal invariance breaking 
in the IHQCD model. 
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1 Introduction 



Heavy ion collision data from the RHIC [T] and LHC |2| experiments suggests that quark-gluon 
plasma near T c behaves as an almost ideal liquid with a very low shear viscosity to entropy ratio, 
rj/s < 0.2. Understanding this small value has turned out to be very difficult using standard 
field theory methods. Calculations in perturbative QCD suggest that this ratio should be 
rj/s ~ 1 |3|, in direct contradiction with the experiments. Moreover, the shear viscosity is a 
real-time transport coefficient, so it is very difficult to calculate in lattice QCD. 



One of the most famous predictions resulting from the AdS/CFT conjecture has been the exis- 
tence of a universal bound for the viscosity to entropy ratio rj/s > l/(Air) |4|, which no physical 
liquid is supposed to violate. Since then, AdS/CFT methods have played an important role 
in increasing our understanding of strongly interacting systems, such as the matter produced 
in heavy ion collisions |5|. However, AdS/CFT is originally a duality between string theory 
and conformally invariant supersymmetric Yang-Mills theory. As QCD is not conformal nor 
supersymmetric, it would be important to develop a holographic model with these features, 
which might be able to describe real- world QCD more closely. One such model is the Improved 
Holographic QCD (IHQCD) proposed by Elias Kiritsis, et. al. j6j [7j, where conformal invari- 
ance is broken by the introduction of a nontrivial potential for the dilaton field, mimicking 
perturbative QCD in the UV region and giving the model a linear glueball spectrum in the IR, 
i.e. making it confining. 



In this talk, we present results published in the papers |8j [9], where using IHQCD we have 
calculated the field theory correlators 

Gf(w,k = 0) = -ijd 4 xe^e(t)([T 12 (t,x),T 12 (0,0)]) (1) 

Gf(w,k = 0) = -i J^xe^%t){[^T^ x )^-T J3 ^M)^ (2) 

where T\ 2 and Tu are the shear and bulk components of the energy-momentum tensor, respec- 
tively. Of our primary interest are the spectral densities, which are defined as 

p s , t ( W ,T)=ImG>,k = 0). (3) 
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Using the Kubo formulas, we can furthermore formulate the relations 

= i im ^ r ) , (4 ) 

C = lira Pb ^ T) (5) 

w^O co 

between the spectral functions and the corresponding transport coefficients, the shear (n) and 
bulk ((■) viscosities of the theory. 

As IHQCD is a two- derivative model, we expect the shear viscosity to entropy ratio to have 
the universal value rj/s = 1/(4-71"), as proven in general in |10] . This is in fact one of our 
reasons to focus on the properties of the spectral densities, as they allow us to distinguish 
between strongly coupled M = 4 SYM and the field theory dual to IHQCD. In the shear 
channel, we furthermore have available a recent pertubative result for the Yang-Mills spectral 
density |11|, which naturally is not applicable in the region of low frequencies, but provides 
us with an important consistency check at higher frequencies. In the bulk channel, the SYM 
spectral density (and bulk viscosity) on the other hand vanish due to conformal invariance, so 
a comparison of our IHQCD results with the perturbative result of |12) and the lattice data of 
|j~3] (for Euclidean imaginary time correlators) is of utmost interest. 



2 Improved Holographic QCD 



Improved Holographic QCD is a five-dimensional bottom-up holographic model, where on the 
gravitational side one starts with a gravity+dilaton system and introduces a potential for the 
dilaton field to describe some crucial features of conformal symmetry breaking |6j[7j. The action 
reads 



S 



16ttGh 



d 5 Xy 



-g 



R 



T 



+ v(\) 



(6) 



where we have used A = e^, with denoting a dilaton field. The form of the potential V(X) is 
found by matching the holographic beta function to the 2-loop perturbative one (in large- N c 
Yang-Mills theory) and by requiring the model to possess a linear glueball spectrum. This can 



be achieved, for example, through the potential 14 



V(X) 



12 
£2 



88 4619 
27 A + W 



+ M1+A) 
1 + A) 2 /3 



(7) 
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Note that this is the potential we have used in our most recent publication |9|, while in |8j a 
slightly different choice, consistent with the original potential of |6|[7), was made. For possible 



issues in the derivation of the potential we would like to point on 15, 16 



A background metric consistent with the above potential reads 

dz 2 



ds z = b z (z) 



-f(z)dt 2 + cix 2 + 



/(*) 



(8) 



where the functions f(z), b(z) and A(z) are determined from the Einstein equations ^ 



W 



1 



1 



AbW 2 - -z{Wf+-bV) 



f 



b = -b 2 W 
3 



A 



f = ZfbW 



(9) 
(10) 

(11) 
(12) 



An additional constraint is that the space-time be asymptotically AdS, i.e. that close to the 
boundary 



6(*)-> 



z ->■ 0. 



(13) 



where C is the curvature radius of the AdS space-time, which without loss of generality we can 



set to C = 1. Using eqs. (|9h-( 12 ), we can next determine the thermodynamic properties of the 



system 17, 18 . Matching these results to the corresponding quantities in large- iV c Yang-Mills 



theory, we are able to finally fix the value of the gravitational constant 



£ 3 



4Ag 
45vr 2 



(14) 



To determine correlators of the energy-momentum tensor on the field theory side, we follow 
(l9l[20| and introduce perturbations around the background metric ((81, 



g 00 = b 2 f (1 + eHoo) , g n = b 2 (1 + eH n ) , g 55 = j (1 + eH 55 ) , g 12 = eb 2 H 12 . (15) 

Here, e is a power counting parameter, and the fluctuation H\2 corresponds to the T\ 2 operator, 
while the fluctuation H\\ is dual to \Tn. Expanding next the Einstein equations to the first 
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order in e, we obtain fluctuation equations for the perturbations 



H 12 + £ \og(b 3 f)H 12 + j 2 H X2 = , (16) 

Hn + £ log(b 3 fX*)H u + - H n = , (17) 

where 

X(A) s (18) 

These equations are to be solved using infalling boundary conditions at the horizon |9|. 

Having obtained the solutions to the above fluctuation equations, the spectral functions are 
available using standard results [19] , 

= fWW (19) 

P ^T) = 6X|fl |gn( ^ 0)|2 . (20) 



In the latter bulk channel result, denotes the quantity of eq. (18) evaluated at the horizon. 
When comparing the holographic results to those of lattice QCD, another useful observable 
turns out to be the directly measurable Euclidean imaginary time correlator. This quantity is 
obtained from the spectral density through the relation 



G(T ,r)=/%( w ,r)^M^lM, j s i/r. (21) 

7T smh (|u;j 







3 Results 



Below, we present our results for the spectral densities, covering the cases of the shear and bulk 
channels separately. 



3.1 Shear Channel 

The ase of the shear channel was considered already in |8j, which, unlike (9], however did not 
make use of the potential of eq. (|7|). Thus, we will here concentrate on the results for the latter 
reference, showing in fig. [T] (left) the ratio of the spectral density and u for various temperatures. 
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Figure 1: Left: the ratio of the IHQCD spectral density p s (oj) and frequency u in the shear channel, 
displayed for three different temperatures in units of C? /(AirGs). The dashed curve cor- 
responds to the asymptotic limit of 7ra; 3 /32. Right: the ratio of the spectral density and 
frequency normalized by the entropy. Using the Kubo formula of eq. we observe that 
the shear viscosity over entropy obtains the universal value of 



We observe a behavior similar to that displayed in fig. 3 of pi; however, due to the choice of 
potential the logarithmic convergence of our result towards its large-w limit (p as = 7rw 4 /32 in 
units where £ 3 /(47rG 5 ) has been set to unity) is somewhat slower. Analogous calculations in 
J\f = 4 SYM theory (cf. [21]) show that for temperatures close to the critical one, a significant 
difference between the IHQCD and SYM results arises, which however rapidly disappears at 
higher temperatures. 



In fig. [T] (right), we next plot p s /u normalized by the entropy, from where we can read off the 
value of the shear viscosity over entropy ratio can as the intercept of the curves at u — > 0. 
As a consequence of IHQCD being a two-derivative model, we find the universal prediction of 
T]/s — l/(47r) to hold at all temperatures considered |10|. 



Despite the limitations of perturbative QCD in the region of low frequencies, at high oj the 
physical behavior of spectral densities is expected to reduce to its predictions due to asymptotic 



freedom 11, 22 [23]. To this end, in fig. [2] (left), we plot both the perturbative (dashed red 
curve) and IHQCD (black curve) results for the shear spectral density, finding that in the UV 
limit, the IHQCD spectral function has the same functional behaviour as the perturbative one 
up to some overall normalization constant. As explained in detail in [9j, this proportionality 
factor can be identified as 4/9 by analytically determining the asymptotic limits of both the 
perturbative and holographic spectral densities. 
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Figure 2: Spectral densities in the shear (left) and bulk (right) channels for T = 3T C , normalized by 
Ad a = A(N^ — 1). The black curves correspond to the IHQCD and the dashed dashed curves 
to the 2-loop perturbative QCD results 11 12 . 



3.2 Bulk Channel 



Next, we move on to present the results of our recent correlator calculations in the bulk channel 
(9). In fig. [3] (left), we first plot pb(cu)/tu for different temperatures, normalized by the factor 
Ad a- Using the Kubo formula (J5|, the bulk viscosity can again be read off as the value of 
the curves in the limit u — y 0. From these results, we see that the bulk viscosity decreases 
with increasing temperature, ultimately approaching the vanishing conformal limit. In fig. [2] 
(right), we next find that in the region of large u, the IHQCD spectral function reduces to the 
perturbative QCD prediction. This reduction to is even more apparent when we study ratio the 
Pb/w 4 plotted in fig. [3] (right), where we see that this ratio quickly approaches the temperature 
independent limit oc 1 / (log co/T c ) 2 . This form in fact coincides with the cu- dependence of the 
field theory gauge coupling g A , to which the leading order perturbative result is proportional 
to 



12, 24 



Finally, lattice simulations of course play an important role in the study of strongly interact- 
ing phenomena, as they are the only fundamentally nonperturbative first principles method 
available. However, due to technical reasons only Euclidean correlators are possible to mea- 
sure on the lattice, and to this end we next specialize to the case of imaginary time correlators, 



cf. eq. (21 ). With this result, it is possible to calculate Euclidean correlators in the both IHQCD 
and pQCD, and furthermore compare the results with the lattice data of |13|, see fig. |4j We 
find that the holographic results seem to be in better agreement with the lattice data than the 
perturbative ones for all temperatures considered. 
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Figure 3: Left: the ratio of the IHQCD bulk spectral density pb(ui) and cj, normalized by 4d^. Right: 
the ratio Pb/u 4 , plotted as a function of uj/(ttT c ) for multiple temperatures, normalized 
by 4d,A- For large values of u, this ratio reduces to the temperature independent limit 
l/(log W /T c ) 2 . 




Figure 4: The imaginary time correlators of the bulk channel, computed for two different temperatures 
in IHQCD (black curves) and pQCD (red dashed curves) |f 2| and compared with the lattice 



data (blue points) of 



13 



4 Conclusions 



In this talk, we have presented the results of a recent calculation of finite temperature correlation 
functions in both the shear and bulk channels of strongly coupled large- N c Yang-Mills theory. 
After motivating our research and stressing the importance of using non-conformal holographic 
models to describe the physics of QCD, we briefly introduced the fundamental properties of 
the IHQCD model and showed how spectral densities can be obtained in this framework. 

Next, we compared our holographic results with corresponding calculations in strongly coupled 
M = 4 supersymmetric Yang-Mills theory |2l] , perturbative Yang-Mills theory 
as lattice simulations in the same theory (13|. In the shear channel, we observed a significant 



11, 121 as well 
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difference between the IHQCD and SYM results close to the critical temperature T c , which 
however rapidly disappears at higher temperatures. We also compared the IHQCD and pertur- 
bative spectral densities in the large-o; region, finding in both cases an asymptotic u 4 behavior, 
albeit with differing overall normalizations. 

Similar methods were next used to calculate the IHQCD spectral density in the bulk channel, 
where we confirmed that with increasing temperature, the bulk viscosity approaches the van- 
ishing result known to occur in conformal theories. Subsequently, we compared our holographic 
calculation with results from perturbative Yang-Mills theory, finding impressive agreement in 
the UV (large-o>)limt, including a perfectly matching overall normalization. We find this ob- 
servation somewhat puzzling, considering the qualitatively different conclusion reached in the 
shear channel. 

Finally, in the bulk channel, we calculated imaginary time correlators using both the IHQCD 
and perturbative QCD spectral functions, and compared the results with available lattice data. 
We conclude that the lattice data seems to favour our holographic results over perturbative one 
for all temperatures considered. 
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